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Abstract 

We present a novel method to search for a new CP-violating phase in the b — > sssss 
transition using B — > (j>(j)X s decays, where X s represents a final state with a specific 
strange flavor such as K , X* 1 * 1 or K^ir^ . Direct CP violation can be enhanced due 
to an interference between an amplitude beyond the standard model (SM) and the 
SM decay amplitude through the r\ c resonance. We find that the CP asymmetry can 
be as large as 0.4 within the present experimental bounds on the b — > sss transition. 
These decays provide a very clean experimental signature and the background is 
expected to be small in particular at e + e~ B factories. A simulation study for 
the B^ — * (fx^K^ decay shows that the statistical significance of CP violation can 
exceed 5 standard deviations with 10 9 B mesons. 
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The phenomenon of CP violation is one of the major unresolved issues in our 
understanding of elementary particles today. In the standard model (SM), CP 
violation arises from an irreducible complex phase, the Kobayashi-Maskawa 
(KM) phase [1], in the Cabibbo-Kobayashi-Maskawa (CKM) weak-interaction 
quark mixing matrix [1,2]. Recent observations of direct CP violation in the 
neutral K meson system [3,4] and mixing-induced CP violation in the neutral 
B meson system [5,6] strongly support the KM mechanism. Despite this suc- 
cess, additional CP- violating phases are inevitable in most of theories involv- 
ing new physics (NP) beyond the SM [7] . Some of them allow large deviations 
from the SM predictions in B meson decays. Examples of such theories include 
super symmetric grand-unified theories with the see-saw mechanism that can 
accommodate the large neutrino mixing [8]. Therefore it is of fundamental 
importance to measure CP asymmetries that are sensitive to the difference 
between the SM and NP. Additional sources of CP violation are also highly 
desirable to understand the origin of the matter-antimatter asymmetry of the 
universe; detailed studies have found no way that CP violation in the SM 
alone could explain the baryogenesis [9]. Many methods to search for a new 
source of CP violation in B meson decays have been proposed until now. One 
of the most promising ways among them is to compare the mixing-induced CP 
asymmetries in the B — > <pKg decay [10], which is dominated by the b — > sss 
transition that is known to be sensitive to possible NP effects, with those in the 
B° — > J/ijjKg decay [11]. Recent measurements by Belle [12] and BaBar [13] 
collaborations yield values smaller than the SM expectation; a difference by 
2.6 standard deviations is obtained when two results are combined. The other 
charmless decays B° — > r/'Kg and B° — > K + K~K° S) which are mediated by 
b — > sss, suu and sdd transitions, also provide additional information [12,13]. 
The present world average with B° — > (f)Kg, rj'Kg and K + K~Kg combined is 
different from the average with B° — > J/ipKg and related modes by 3.1 stan- 
dard deviations [14] . Possible implications of these measurements are discussed 
in the literature [15]. 

In this Letter, we propose a novel method to search for a new CP- violating 
phase in the hadronic b — > s transition using B ± — > (fxfiXf decays [16]. Here 
Xf represents a final state with a specific strange flavor such as K ± or K^. 
These non-resonant direct decay amplitudes are dominated by the b — > sssss 
transition. A contribution from the b — > uus transition followed by rescattering 
into sss is expected to be small because of the CKM suppression and the 
OZI rule [17]. In these decays, when the invariant mass of the (fxf) system is 
within the r] c resonance region, they interfere with the B^ — > > tpcfyXf 
decay that is dominated by the b — > ccs transition. The decay width of r\ c 
is sufficiently large [18,19] to provide a sizable interference. Within the SM, 
this interference does not cause sizable direct CP violation because there is 
no weak phase difference between the b — + sssss and the b — > ccs transitions. 
On the other hand, a NP contribution with a new CP- violating phase can 
create a large weak phase difference. Thus large CP asymmetries can appear 
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only from NP amplitudes, and an observation of direct CP violation in these 
decays is an unambiguous manifestation of physics beyond the SM. 

Although the same argument so far is applicable to the B ± — > 4>Xf decays, 
there is no guaranteed strong phase difference that is calculable reliably for 
these decays. In contrast, the Breit-Wigner resonance provides the maximal 
strong phase difference in the case of B ± — > ((f)(f)) m ^, mr]c Xf decays. CP vi- 
olation that arises from interference between a resonant and a nonresonant 
amplitude was initially studied in top-quark decays [20] and in radiative B ± 
decays [21]. The charmonium width effect on the SM direct CP asymmetries 
in B decays was first considered in 5 ± — > ^(Xco) 71 " 1 * 1 decays [22] and was also 
used in subsequent theoretical studies [23]. 

Experimentally, an evidence for the decay B^ — > r\ c K ± was first reported by 
the CLEO collaboration [24]. Recently Belle [19] has also reported the first 
observation of the B° — > i] c K*° decay. This implies that other modes such as 
B + — > i] c K* + will also be seen with a similar branching fraction, so that we 
will be able to study semi-inclusive B ± — > r) c Xf transitions experimentally. 
The semi-inclusive branching fraction of B ± — > r] c Xf is not yet measured, but 
is theoretically expected to be comparable to the branching fraction of the 
semi-inclusive decay B ± — > J/ipXf [25]. 

We derive the rates and the asymmetry of the decays B^ — >■ (4>4>) m ^ mric Xf 
based on the formalism described in the study of B* 1 — * rj^Xco)^ decays [22]. 
The distribution of two 0's is determined with two kinematical variables; one 
is the invariant mass of the 00 system, to, and the other is the angle 9 between 
the 5-meson momentum and the momentum of one of two 0's in the center-of- 
mass frame of the 00 system. To have the interference between resonant and 
direct amplitudes, to should be in the r\ c resonance region. To be specific, we 
require in this study that the difference between to and 7] c mass (M) should 
satisfy \m — M\ < 3r, where V is the the width of the r] c resonance [26]. The 
differential decay rate normalized with the total B^ decay rate is then given 
by the following equation: 



where R(s) is the resonant amplitude, -D ± (s, z) is the direct amplitude of the 
B^ — > (jxpXf decay, T B is the total B decay rate, s = to 2 and z = cos 9. 

The resonant amplitude R(s) is given by 
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where or is a product of the weak decay amplitude of B — > rj c Xf and the 
real part of the 77 c decay amplitude to 00. 

The direct amplitude .D 1 * 1 is separated into contributions from the SM, D S m, 
and from NP, D^ P , 



D ± (s « M 2 , *) = D S m(s « M 2 , z) 

+ D± P (s^M 2 ,z), (3) 



D S M( S -M 2 ,,) = ^e^ (4) 
D&(8*ltf,z) = ^e«e»°«*, (5) 

where od(z) is a real part of the SM direct amplitude, 8 (8') is a strong phase 
difference between the resonant amplitude and the SM (NP) direct amplitude, 
atfp(z) is a real part of the NP amplitude and Onp is a new CP- violating 
phase. If 8 ^ 8' holds, direct CP violation can also occur from an interference 
between the SM and NP direct amplitudes. We do not take this case in our 
study and assume 8 — 8' in the following discussion. 

The difference between the decay rates of B + and B~ is given by 

— (— — ) =7~(z) = -4ira R a NP (z) coso" • sin6 N p . (6) 

1 r dz dz 



Similarly the sum of two decay rates is given by 



i c/r^~ (/r — 

_(__ + __) = 7 +( z ) ^ 2nal + 24a 2 , (z)(r 2 + 2r cos 6 N p + 1) 

— 47raRaD(<2)(V cos ©np + 1) sin 5, (7) 

where r = aNp(z)/an(z) is the amplitude ratio of NP to the SM. The z 
dependence of r reflects the spin components of the 00 system, which can be 
determined at B factories in the future from the differential decay rates in 
the mass-sideband region below the r] c resonance. Although only a pseudo- 
sealer component in the direct transition interferes with the r] c resonance, the 
effect of other components can be estimated by such a measurement. Thus, 
for simplicity, we assume that the direct transition is dominated by a pseudo- 
sealer component and ignore the z dependence of r in the following discussion. 

The maximum asymmetry is realized when cos 8 ~ 1 is satisfied. Assuming 
that 8 is small following the discussion by Eilam, Gronau and Mendel [22], 
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the differential partial rate asymmetry is 



A ( ) = 7 (jj ^ -47ra R g N p(z) sin 6 NP 

CPlZj ~ 7+(*) ~ 2vra| + 24a^(z)(r 2 + 2r cosGnp + 1) ' 1 J 



As a measure of CP violation, we define the following CP-asymmetry param- 
eter: 



Acp = 



f\dzj (z) 2 



The numerator of Aqp can be expressed with the branching fraction of the 
resonance (27rapJ and that of NP in the resonance region (£>np): 

i 

r 2tt 
J d Z1 -{zf = (2vr4) • B NP • — sin 2 6 NP , (10) 

-l 



where 

2tt4 = B(B ± - n c Xf) ■ B( Vc - 00) , (11) 

and 

(M+3r) 2 i 

Bnp =mt J ds j dza lM- ( 12 ) 
(M-3r) 2 -i 

The inclusive branching fraction B(B ± — > r] c Xf) has not been measured and 
only an upper limit of 0.9% is available [18]. Therefore we take theoretical 
expectations that range from 0.3% to 0.7%, depending on the estimation of 
the ratio of B(B -> r? c X) to B(B -> J/^A") [25]. For £(?? c -> 00), we take the 
present world average of (0.71 ± 0.28)% [18]. As a result we obtain 2jca 2 l = 
(2~5)xl0- 5 . 

There is no direct experimental bound on the B — > <f)(f)X s decay (see Note 
addedin the end of this Letter). A bound on the B — > (f>X s decay (2.2 x 10~ 4 at 
90%C.L.) is available [27]. The search, however, was devoted only for energetic 
mesons above > 2.1 GeV/c. On the other hand, the momentum of 
mesons in the B — > <f)<f)X s decay is lower; it should be less than 2.2 GeV/c. 
Consideration on other experimental bounds also leads to a weak limitation 
on b — > sg [28]; B(b — > sss) ~ 1% is still allowed, while the estimation of 
B(b — > sss) within the SM is ~ 0.2% [29]. Since there is no reliable way to 
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calculate the hadronic matrix element of the multi-body decay B — > 4>(f)X s , we 
use an event generator that is based on the LUND fragmentation model [30] for 
the estimation of <Bnp- We start from the b — > sg* — > sss transition that mainly 
produces soft mesons. We find that it is possible to assume B(b — > sg* — > 
sss) ~ 1% without conflicting the existing experimental results such as the 
branching ratio of the exclusive B — > decay [31]. We use the default value 
of the ss popping probability in JETSET7.4 [30], which is consistent with the 
ratio of the branching fractions between B — > J/ipK<j) and B — > J/ipK [32]. 
We estimate £> NP to be ~ 5 x 10 -6 for Bib — > sg* — > sss) = 1%. Estimations 
of the corresponding branching fraction within the SM typically yield smaller 
values than this estimation. For example, the branching fraction of B — > <f)X s 
is expected to be around 10 -4 . With an additional ss popping and with the 
aforementioned simulation, we obtain B(B — > (p(f>X s )s M ~ 3 x 1CT 7 in the 
i] c mass region. Also from the expected B(b — > sss) in the SM, we obtain 
~ 9 x 10~ 7 . Thus r 2 ~ 5 holds for B(b -> s#* -> sss) = 1%. Note that, even 
in this case, an observation of non-resonant B — > (fxf)X s decays in the i] c mass 
sideband region alone can not establish physics beyond the SM because of the 
difficulty in estimating branching fractions of multi-body hadronic decays. 

The evaluation of the denominator of Aqp is straightforward with the esti- 
mations mentioned above. We obtain the following result for Bib — > sg* — > 
sss) ~ 1% and 2na 2 K = 2 x 1(T 5 : 



A large CP asymmetry of 0.4 is allowed. The asymmetry is roughly propor- 
tional to |r|. Therefore it can be sizable even with r 2 < 1; for example, Acp ~ 
0.1 is allowed for r 2 = 0.3. 

A CP asymmetry proportional to sin0 3 (<f> 3 = 59° ± 13° [18]) may arise 
within the SM due to an interference between the resonant amplitude and a 
SM contribution from the b — > uus transition followed by rescattering into 
sss. The effect is expected to be small since the b — > uus decay ampli- 
tude is CKM-suppressed, and the rescattering amplitude is also suppressed 
by the OZI rule. From the measured branching fractions for B + — > K + K~ K + 
and B + — > K + K~n + together with the additional O(A) suppression, where 
A = 0.2229 ± 0.0022 is the sine of the Cabibbo angle [18], the fraction of the 
b — > uus transition in the B + — > K + K~ K + decay amplitude is estimated to be 
V0.022 ± 0.005 in Ref. [33]. Using this information and with a modest assump- 
tion on a rescattering amplitude ratio Aiuu — > (p(p)/A(uu — > K + K~) < 0.3 
from the OZI rule, we estimate the CP asymmetry to be less than 0.9%. Thus 
we conclude that the CP asymmetry within the SM is most likely below 1%. 



A CP = 
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We perform Monte Carlo simulation for the B — > (fxf)K decay and estimate 
statistical errors on the CP asymmetry parameter. For this decay mode, the 
background level is expected to be small enough to be neglected [34]. The 
reconstruction efficiency and the 00 mass resolution are estimated using a 
GEANT-based detector simulator for the Belle detector [35]. We obtain ~300 
events for Ng = 10 9 , where Ng is the number of charged B mesons recorded by 
a detector. We perform an unbinned maximum-likelihood fit to the differential 
decay rate distribution, which is proportional to \R(s) + D ± (s, z)\ 2 , instead of 
integrating the distribution. We choose the following two free parameters in 
the fit: Aq P = — 2r(a D /a R ) sinO NP and£> = a^(r 2 +2r cos NP +1). A° CP is the 
CP asymmetry in the Breit-Wigner term. B is proportional to the branching 
ratio of the non-resonant B ± — > (jxfiK^ decay below the i] c mass region. The 
statistical error for A Q CP is estimated to be 5A° CP ~ 0.06. Figure. 1 shows the 
5a search regions for Ng = 10 9 (dotted line) and for N B = 10 10 (solid line), 
which will be accessible at next-generation high-luminosity e + e~ B factories. 
Direct CP violation will be observed in a large parameter space above a 5er 
significance. 

The new CP- violating phase 9np also affects time-dependent CP- violating 
asymmetries A C p(t) = S sm(Am d t) + Acos(Am d t) in B° — > 4>K° S and related 
decays. Here Am d is the mass difference between the two B° mass eigenstates, 
and S and A are parameters for mixing-induced CP violation and direct 
CP violation, respectively. Ignoring a strong phase difference between the 
amplitude of NP (A NP ) and SM (Asm), we obtain 

sin 20x + 2psin(20! + 9 NP ) + p 2 sin(20! + 29 NP ) 
~ l + p 2 + 2pcos e NP ' ^ ' 



where p = A^p/Asm is an amplitude ratio of NP to the SM and 0i is one of 
the angles of the unitarity triangle. In particular, a difference in S between 
B° -> <j>K% and B° -> J/ipK° s decays, i.e. AS = S(<pK° s ) - S(J/ipK° s ) ^ 0, 
would be a clear signal of the new phase since S(J/ipKg) = sin20i is held 
to a good approximation. We define expected statistical significance of the 
deviation from the SM by A° CP / 5A° CP for the B^ 1 — > $$K ± decay and by 
AS I b AS for the B° — > 4>K° S decay, where 5 AS is an expected statistical 
error of AS extrapolated from the latest result by the Belle experiment [12]. 
Although r 2 is not necessarily equal to p 2 , both decays are governed by the 
same b — > sss transition. Therefore it is reasonable to choose r 2 = p 2 for 
comparison. Figure. 2 shows the resulting significance for 10 10 B mesons and 
with r 2 = p 2 = 0.5. The significance for AS largely depends on the sign of 
6 NP , which is not the case for the B ± — > 00_ft' ± decay. The sign dependence 
arises from an asymmetric range for AS; to a good approximation, we have 
-l-sin20i < AS < 1 -sin 20i where sin20i = +0.736±0.049 [14]. Therefore 
the B ± — + (fxpK^ 1 decay plays a unique role in searching for a new CP- violating 
phase. 
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In the above estimation, we use parameters that have uncertainties. However, 
they can in principle be measured precisely if a sufficient number of B mesons 
are produced. In our estimation, we assume efficiencies and background lev- 
els that have been achieved with the Belle detector at the KEK B factory. 
They depend on the actual detector performance and beam conditions, which 
might be different at a next-generation B factory with the higher luminosity. 
Detailed simulation studies as well as some extrapolation from data at current 
B factories will be needed for further quantitative evaluation. 

Experimental sensitivities can be improved by adding more final states. The 
technique to reconstruct X s , which has been successfully adopted for the mea- 
surements of semi-inclusive B — > X s ££ transitions [36] , can be used for this pur- 
pose. Flavor-specific neutral B meson decays, such as B° — > 00fT*°(^ K + n~), 
and other charmonia such as the Xco — > 00 decay can also be included. The 
method proposed in this Letter is also applicable to other beauty hadrons. 
Examples include B s — > rj c {-^ 00)0 an d Aft — > rj c {-^ 00) A. 

In summary, we have investigated B — > 00X S decays as a probe for a new 
CP- violating phase in the b — > sssss transition. Large CP violation may arise 
from an interference between B + — > i] c (^ 00) and a contribution of the 
physics beyond the SM to the non- resonant decay B + — > 00X+. We find that 
the direct CP asymmetry can be as large as 0.4. The experimental signature 
will be very clean in particular at e + e" B factories; when 10 9 charged B 
mesons are available, the statistical error on the CP asymmetry parameter is 
estimated to be ~0.06. The statistical significance of CP violation in the B^ —± 
<fxf)K ± decay can exceed 5 standard deviations. An even better sensitivity will 
be expected by including other final states. We also find that this result is 
comparable to the expected sensitivity with mixing-induced CP violation in 
the B° — > 4>K° S decay. Thus B — > 00X S decays will be useful to search for CP 
violation beyond the SM in the b — > sssss transition at high-luminosity e + e~ 
B factories. 
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Note added. -After finishing this study and as we were preparing the manuscript, 
the Belle Collaboration announced evidence for B — > cjxpK decays [37]. The 
signal purity is close to 100% when the 00 invariant mass is within the r\ c 
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mass region, as assumed in our study. The results are still limited by small 
statistics; using the measured branching fractions of the B — > <pcf)K decay 
within and below the r\ c mass region, we find that the maximal CP asymme- 
try of Acp ~ 0.4 is still allowed. The same paper also reports B{r\ c — > 00) = 
(1.8io;g =t 0-7) x 1CT 3 , which is smaller than the current world average. We 
repeat the fit procedure described in this Letter with the above branching frac- 
tion. Although the smaller value for B(rj c — > 00) results in the smaller number 
of signal events, the CP asymmetry from the interference between the reso- 
nant and the new-physics amplitudes becomes larger. We find that the change 
in B(r] c — ► 00) does not largely affect the significance; the difference is less 
than 10% for r 2 = 0.5 and suiOnp = 1. 
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Fig. 1. Expected sensitivities on direct CP violation in the — ► (fxfrK^ decay for 
10 9 B mesons (dotted line) and 10 10 B mesons (solid line). In the regions above the 
curves, direct CP violation can be measured with a 5a significance or larger. 




NP (degree) 

Fig. 2. Expected statistical significance of deviations from the SM for direct CP 
violation in the B^ — > cfxpK^ decay with r 2 = 0.5 (solid line) and for time-dependent 
CP violation in the B° — ► <pKg decay with |t1np/^4sm| 2 = 0.5 (dashed line). For 
each case, significance is calculated with 10 10 B mesons. 
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